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ABSTRACT 

We present an observational study of the planetary nebula (PN) NGC 7354 consisting of 
narrow band Ha and [N II]A6584 imaging as well as low and high dispersion long-slit spectroscopy 
and VLA-D radio continuum. According to our imaging and spectroscopic data, NGC 7354 has 
four main structures: a quite round outer shell and an elliptical inner shell, a collection of low- 
excitation bright knots roughly concentrated on the equatorial region of the nebula and two 
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asymmetrical jet-like features, not aligned neither with the shells axes, nor with each other. 
We have obtained physical parameters like electron temperature and electron density as well as 
ionic and elemental abundances for these different structures. Electron temperature and electron 
density slightly vary throughout the nebula going from ~ 11, 000 to ~ 14, 000 K, and from ~ 1000 
to ~ 3000 cm~'^, respectively. The local extinction coefhcient cjj^ shows an increasing gradient 
from South to North and a decreasing gradient from East to West consistent with the number 
of equatorial bright knots present in each direction. Abundance values show slight internal 
variations but most of them are within the estimated uncertainties. In general, abundance values 
are in good agreement with the ones expected for PNe. Radio continuum data are consistent 
with optically thin thermal emission. Mean physical parameters derived from the radio emission 
are electron density Ue = 710 cm"^ and M(HII)= 0.22 Mq. 

We have used the interactive three-dimensional modeling tool shape to reproduce the ob- 
served morphokinematic structures in NGC 7354 with different geometrical components. Our 
observations and model show evidence that the outer shell is moving faster (~ 35 kms~^) than 
the inner one ~ 30 kms~^. Our SHAPE model includes several small spheres placed on the outer 
shell wall to reproduce to equatorial bright knots. Observed and modeled velocity for these 
spheres lies between the inner and outer shells velocity values. The two jet-like features were 
modeled as two thin cylinders moving at a radial velocity of ~60kms^^. In general, our shape 
model is in very good agreement with our imaging and spectroscopic observations. Finally, after 
modeling NGC 7354 with shape, we suggest a possible scenario for the formation of the nebula. 

Subject headings: planetary nebula: individual: NGC 7354 — ISM: kinematics and dynamics — ISM: 
abundances 



1. Introduction 

Since the recognition that winds from the 
central stars of planetary nebulae (PNe) play 
an import ant role in the shaping of thes e ob- 
jects fe.g. iKwok et al. Ill978t iBahckl [l987h . and 
that only a small fraction of them show circu- 
lar symmetry, the interest in PN morphology 
has triggered an active field in both theoreti- 
cal and observational astronomy. Many studies 
have been devoted to classify (e.g., Balickl 1987 ; 
Schwarz et al. Ill992l: iManchado et al. Ill996[) a^ 



to model the basic morphologies observed (cir- 
cular, elliptic al and bipo l ar) w i th noticeable 



succe ss (e.g., iBalick et al.l 119871 : iHajian et al. 



19971 ). However, high resolution images have 
shown that the morphologies of PNe are far from 
simple. Multiple shells, multipolar structures, 
highly coUimated outflows, microstructures and 
peculiar geometries are present in many PNe 
and cani iot be explained with simplistic mod- 
els (e.g., 'Sahai & Trauger "l998"; 'Mira nda et"ar 
[2OO6 : Miranda, Pcrcira & Guerrero 2009|). Nowa 
days, it is accepted that highly coUimated out- 
flows play a crucial role in the shaping of PNe 
( Sahai fc Trauger 1998h . Nevertheless, other 



physical processes are probably present, so that 



most likely the shaping of PNe is a result of 
many processes actin g at the same time (e.g., 
Bahck fc Fran^ |2002| ). In order to understand 
complex PNe, the first step is to identify the 
structural components present and to define their 
nature. High-resolution, spatially resolved spec- 
troscopy combined with narrow-band imaging has 
demonstrated to be a powerful tool to disentagle 
the structural components in PNe, to infer their 
nature and to constraint t he ejection processes in- 



volved in their formation ( Miranda fc Solf 



Vazauez et al. 


1999 


VazQuez et al. 


1998 



2008 



Guerrero et al. 



Lopez et al. 1120001 ). 



1992; 



2008 : 



Although at first glance NGC 7354 looks like 
an elliptical PN, extense Ha, [O III]A5007, 
[N II]A6584 and mid-infrared imaging and spec- 
troscopic data have revealed that it p ossess a more 



comp l ex structure (lSabba,din et al. 



19871 : iHaiian et al. 



997 



1983; Balic 
Phillips erai. I f200 



In particular, iBalickl (jl987[ ) described this nebula 
as consisting of an inner halo, a thin bright rim, 
two spike-like tails, and low-excitation patches 
projected onto the rim-halo interface. The large- 
scale structures observed in NGC 7354 are quali- 
tatively well described using the i nteracting winds 



theory (e.g. iBalick et all Il987t iMellema I Il995l: 
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Haiian et al. 1120071 ) but no deep analysis has been 
carried out for the smah-scale structures and their 
relationship with the large-scale ones. 

In this work we present a detailed observational 
study of the morphology, internal kinematics, and 
physical and chemical properties of NGC 7354. 
These data allow us to discuss and model each of 
the components present in the object and to sug- 
gest a possible scenario for the formation of this 
nebula. 

2. Observations and Results 
2.1. Optical Imaging 

Narrow-band images of NGC 7354 were ob- 
tained on 1997 July 24 with the Nordic Opti- 
cal Telescope (NOTfl and the HiRAC camera 
equiped with a Loral CDD of 2048x2048 pixels 
and a plate scale of 0"11 pixel^^. Two narrow- 
band fihers were used: [N II]A6584 (AA = 10 A) 
and Ha (AA = 10 A). The exposure time for both 
filters was 900s. Seeing was about 0"9 during the 
observations. 

Figure [1] shows a mosaic of our Ha and [N II] 
images, including unsharp masking images in both 
filters constructed to show up both the large- and 
small-scale structures in the nebula. In these im- 
ages, we can identify the main structures pre- 
viously described by other authors: the outer 
shell, the elliptical inner shell, the bright equa- 
torial knots, and the two jet-like features. In the 
following we will describe in more detail each of 
these structures as well as new morphological de- 
tails that have not been previously mentioned. 

The outer shell looks like a round envelope in 
the high-contrast images but it appears as a faint 
cylindrical structure in the low-contrast images. 
Its size is ~ 33" x 29" with the major axis oriented 
at position angle PA~ 15°. The outer shell is 
brighter in Ha than in [N II] although in [N II] 
several bright knots are observed at the edges of 
it. 

The inner shell present an elliptical shape with 
its major axis oriented at PA ~ 30° and a major 
and minor axis length of ~ 21' x 16", respectively. 



This structure is noticeable fainter in [N II] than 
in Ha. Moreover, the regions along the minor axis 
are particularly bright in Ha. A closer inspection 
of the images shows that the polar regions deviate 
from the elliptical shape and appear as two bub- 
bles. This is particularly noticeable in Ha. We 
will refer to these regions of the inner shell as the 
polar caps . 

The bright equatorial knots are observed in 
[N II] but not in Ha. They are mainly concen- 
trated in two groups. East and West, and along 
the equatorial plane of the outer and inner shells. 

The two jet-like features are bright in [N II] but 
much fainter in Ha. The northern feature is ori- 
ented at PA~ 13° and extends ~ 7". The southern 
feature is oriented at PA ~ 205°, extends ~ 11", 
and appears narrower than the northern one. We 
note that the orientation of these two features does 
not coincide with each other nor with the orienta- 
tion of the outer and inner shells. 

In order to improve the view of NGC 7354, we 
retrieved a Hubble Space Telescope (HST) image 
from the MAST ArchiveH (Proposal ID: 7501; P.I.: 
A. Hajian; date of observation: 1998 July 21; fil- 
ter F658N; exposure time lOOOsec). Fig. [5]shows 
this image. The morphology of the outer and in- 
ner shells is similar to that observed in ground 
based images (Fig. 1). The jet- like freatures ap- 
pear as cometary tails with a bright knot facing 
the central star and faint tails directed outwards. 
They present a knotty structure, particularly the 
southern one. The HST image resolve the bright 
equatorial knots into a series of small knots and 
filaments embedded in diffuse emission. 

2.2. Radio continuum 

The A3. 6 cm radio continuum observations were 
made with the Very Large Array (VLA) of the 
NRACH on 1996 May 31 in the DnC configura- 
tion. The standard VLA continuum mode with 



^The NOT is operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Institute de Astrofisica de Canarias. 



^Some of the data presented in this paper were obtained 
from the Multimission Archive at the Space Telescope Sci- 
ence Institute (MAST). STScI is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc,, under 
NASA contract NAS5-26555. Support for MAST for non- 
HST data provided by the NASA Office of Space Science 
via grant NAG5-7584 and by other grants and contracts. 

^The National Radio Astronomy Observatory (NRAO) is 
operated by Associated Universities Inc., under cooperative 
agreement with the National Science Foundation. 
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a bandwidth of 100 MHz and two circular polar- 
izations was employed. Phase center was set at 
RA(2000.0)=22M0'"20!1,DEC(2000.0)=+61°17'06'.'0. 
Flux calibrator was 3C48 (adopted flux density 3.3 
Jy) and phase calibrator was 0019+734 (observed 
flux density 1.0 Jy). Total on-source integration 
time was 30 minutes. The data were calibrated 
and processed using standard procedures of the 
Astronomical Image Processing System (AIPS) 
package of the NRAO. 

Fig. [3] shows an uniform- weighted map of 
NGC 7354. The emission presents a circular mor- 
phology and extends ~42" in diameter. Two emis- 
sion maxima are observed separated by ~10"and 
oriented at PA~ 100°. These emission maxima co- 
incide with the Ha bright regions observed along 
the minor axis of the inner shell. 

From our data we derive a peak flux density of 
72mJybeam-i at position RA(2000.0)=22M0'°2ff:44, 
DEC(2000.0)=-h61°17'06'.'8, and a total flux den- 
sity of 502 ± 4 mJy. We note that our map i s 
similar to that obtained bv lTerzian et al. " (1974), 
being the flux density values obtained in both 
works consistent with each other and with op- 
tically thin thermal emission. Adopting a dis- 
tance of 1.5 kpc for the nebula ( Sabbadin 19861 
see Sec. 3.1) and fo ll owing the formulation by 
Mezger fc Henderson I (jl967l ). we derive a mean 
electron density of ~ 710 cm~^ and an ionized 
mass of ~ 0.22 A/q. 

2.3. Low Resolution Spectroscopy 

Low resolution, long-slit spectra were obtained 
with the Boiler & Chivens spectrograph mounted 
on the 2.1m telescope at the San Pedro Martir 
Observatory (0AN-SPMj3 during three observing 
runs: 2002 June 14, 2002 August 7, and 2002 De- 
cember 10 and 11. A CCD SITe with 1024 x 1024 
pixels was used as a detector. We have used a 
400 lines/mm dispersion grating and a slit width 
of 2" giving a spectral resolution (FWHM) of 7A. 
Spectra reduction was performed using 



*The Observatorio Astronomico Nacional at San Pedro 
Martir (OAN-SPM) is operated by the Instituto de 
Astronomfa of the Universidad Nacional Autonoma de 
Mexico. 

^The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement 



standard procedures. We have used four slit posi- 
tions to cover specific regions of the nebula. Fig. U] 
shows these slits, labelled A to D, and the regions 
extracted from each long-slit spectrum overplot- 
ted on the unsharp masking [N II] image. Regions 
are denoted by a letter refering to the slit position 
followed by a sequence number along the slit. In 
total, 21 regions in NGC 7354 have been analyzed. 

Table 1 presents dereddened line fluxes, ob- 
served H/3 flux and the logarithmic extinction 
coeficiente ch^ for each region, the last has 
been derived from the observed Ha/H/3 ratio as- 
suming case B recombination. The fluxes have 
be en dereddened aplying t he extinction law by 
iCardelh. Clavton fc Mathid (fl989i) . Electron tem- 
perature (Te) has been derived from the [O III] 
and/or [N II] emission lines, while electron den- 
sity (iVe) has been derived from the [S II], [CI III] 
and/or [Ar IV] emission lines, in both cases 
using the Five-Lev el Atom Diagnostic Package 
NEBULAR tDeRobertis. Dufour fc Hunt I Il987t 
Shaw fc Dufour Ill994ll from IRAF. Derived phys- 



ical parameters and their estimated errors are 
shown in Table [H Error estimates take into ac- 
count the readout and photon noise and they are 
propagated along the calculation of physical quan- 
tities. Whenever possible, Ng values were calcu- 
lated using the value of Tg derived from the cor- 
responding high- or low-excitation ion, [O III] or 
[N II] , respectively. 

Extinction within the nebula, as indicated 
by Cup, shows a slight increase from South to 
North with values ranging from ~1.7 to ~2.4, 
respectively. Relatively high values of ch/3 ~ 2.4 
are found in the low-excitation equatorial knots. 
These values can be compared to those in the 
surroundings of the knots where the extinction 
decreases up to cnp — 1-7. This result suggests 
that the equatorial knots are denser and/or con- 
tain more dust than the rest of the nebula. 

Electron density slightly increases inwards from 
the outer shell (regions Bl and C6) with values 
of ~ lOOOcm^'^ to the inner regions of the in- 
ner shell (region B4) where values of ~ 2 400 cm~^ 
are found. The bright equatorial knots are clearly 
denser that both the outer and inner shells with 
an average density of ~2 600cm~^. For the jet- 
like features, the electron density is relatively low 



with the National Science Foundation. 
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with values around ~ 1 300 cm 

Electron temperature also seems to slightly in- 
crease from ~ 13,000 K at the walls of the outer 
shell (regions Bl, B6 and B7) to ~ 15,000 K at 
the center of the inner shell (region B3). In 
the bright equatorial knots, electron temperature 
ranges from ~ 10,000 to ~ 12,400 K. For the north- 
ern jet-like feature, no [O III] and [N II] lines were 
detected with a good signal-to-noise; in the south- 
ern jet-like feature, we have obtained an electron 
temperature of ~ 12,000 K (regions A5 and D3). 

Ionic and elemental abundances values are 
listed in Tables 3 and 4, respectively. They have 
been obtained with the task IONIC in IRAF 
and using the ionizat i on co rrection factors by 
Kingsburgh fc Barlow I (|1994l ). respectively. Ta- 
ble m also lists abundances in other objects for 
comparison purposes. Small abundance vari- 
ations, within the estimated uncertainties, are 
observed throughout NGC 7354. In general, we 
found that all our abundance determinations are 
consistent with those of PNe, see Table ID Since 
neither He nor N overabundance is observed, as in 
the case of Type I PNe, we can say that it behaves 
like a Type II PN. 

Elemental abundances for NGC 7354 have been 



reported by severa l authors (IHaiian et al 



Perinotto et al. 



1997 



2004al : 



Martins fc Viegas 200 

Stanghellini et aLll2006( ). However, a comparison 
of our abundance estimations with those from the 
literature is not straightforward since we have ob- 
tained abundance values for several slit positions 
and specific r e gions along them. In the case of 
Hajian et al~l (|l997t ). they report abundance val- 
ues of six regions along a slit position very sim- 
ilar to our slit A. Their abundances tend to be 
two or three times higher than ours in the case of 
0/H, N/H and S/H but lower in He/H and sim- 
ilar in Ar/H. Comparison with other studies can 
only be done through average values along each 
of our slit positions or the total average obtained 
from all our studied regions. We have compared 
our total averaged Q/H N/H, S/H, Ar/H abun- 



dan ces with those froni iPerinotto et al. I (j2004al ) 
and Stanghellini et al."l (|2006[ ). In both studies, 
their abundan ces seem to be lowe r thari ours but 
in the case of IStanghellini et al. ] (|2006l) this dif- 



2.4. High Resolution Spectroscopy 

High-resolution, long-slit spectra were obtained 
in 2002 July 15 to 17 and 2007 July 10 to 17 
with the Manchester Echelle Spectrometer (MES; 
(|Meaburn et al.ll2003l )) mounted on the 2.1m tele- 
scope at San Pedro Martir Observatory (OAN- 
SPM). A Site CCD with 1024x1024 pixels was 
used as a detector. Binnings of 1 x 1 and 2x2 
were used in 2002 and 2007, respectively. Slit 
width was 1.6" and the achieved spectral resolu- 
tion (FWHM) is 12 km s"^ The slit was oriented 
North-South and centered at several right ascen- 
sions across the nebula, except in those slits that 
covered the jet-like features, which have been cen- 
tered on the central star and oriented at PAs 13° 
and 25°. Figure 5 shows the used slit positions, 
numbered from 1 to 8, superimposed on the un- 
sharp masking [N II] image. Data reduction was 
carried out with the IRAF package. 

Grey-scale, position-velocity (PV) maps, de- 
rived from the eight long-slit spectra, are shown 
in Fig. E] for the Ha, [N II]A6584, and He IIA6560 
emission lines. From the long-slit spectra, we de- 
rived a heliocentric systemic velocity of —42 ± 
2kms~^ for NGC 7354, in exc ellent agreement 
with —41 ± 2kms^^ deduced bv lSabbadin et al. 



ference may be due to tha fact that they consider 
the PN as a whole and excluded special features 
in the nebula. 



(|l983l ). Through this paper, we will consider the 
systemic velocity as the origin for quoting internal 
radial velocities and the declination of the central 
star as the origin for quoting distance measure- 
ments. The PV maps allow us to recognize the 
structural components that have been identified 
in the direct images. In the following, we will de- 
scribe the spatio-kinematical properties of these 
components. 

The [N ii] emission from the outer shell is recog- 
nizable at slit positions 1,2,3,5, and 6. The emis- 
sion is very faint and in the central nebular re- 
gions it appears superposed by the stronger emis- 
sion from the inner shell. Maximum radial velocity 
of ~ 35 kms^^ is observed at the center and de- 
creases with distance to the central star. Emission 
from the outer shell can be recognized in the Ha 
line (Fig. 6) by its spatial extend. However, the 
large thermal width and, probably, low expansion 
velocity in this line do not allow us to obtain de- 
tailed kinematic information. The outer shell can- 
not be recognized in the long-slit spectra of the 
He IIA6560 line. 
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The inner shell can be identified at all slit po- 
sitions in the three lines. The emission feature of 
the inner shell appears as a velocity ellipse in the 
PV maps. The size of the velocity ellipse is maxi- 
mum at slit 8 with values of 30" in Ha and [N II] , 
and 20" in He II. Maximum velocity splitting is 
observed at the center of the nebula (slits 3, 7, 
and 8) and amounts 60, 56 and 48 kms~^ in Ha, 
[N II], and He II, respectively. This implies expan- 
sion velocities for the inner shell of 30, 28 and 24 
kms""'^, respectively. Since we estimate an error 
of < Ikms"""^ in our velocity measurements, these 
results show that the He II emission is confined 
to a slow expanding, inner thin layer of the inner 
shell while the [N II] traces the outer layer that 
expands faster. 

It is worth noticing the very faint and wide 
emission observed in [N II] (slit positions 2, 3, 7 
and 8) located at about 12" and 8" to the North 
and South of the central star, respectively. The 
velocity of this faint emission spreads from ~ —45 
to ~ 33kms^^. Comparing our [N II] direct image 
and these particular PV maps, we identify these 
weak emissions with the two polar caps located 
on the main symmetry axis of the inner elliptical 
shell. 

All slit positions in our [N II] PV maps, show 
the presence of at least one of the bright low- 
ionization knots identified in our optical [N II] im- 
ages (Fig. [11 right panel). In all these PV maps we 
can see that most of them are located very close to 
the zero position line, i.e. almost on the equato- 
rial plane of the nebula. Although these PV maps 
show that emission arising from different bright 
knots may be mixed due to projection effects, we 
can estimate that their expansion velocity ranges 
between the inner and outer shell expansion ve- 
locities, 28kms~^ to ~40kms^^, with only two 
of them having lower velocities, ~24kms^^. A 
final remark on the low-ionization knots is that 
in our [N II] PV maps, slit positions 1 and 6, we 
can see two weak spots of emission almost simetri- 
cally located at about 14" to the North and South 
from the central star, respectively. Both emis- 
sion spots show negative and very low velocities 
of ~ — 5kms~^. We identify these two spots of 
emission as coming from two bright knots that, in 
projection, appear to be located on the outer shell 
wall (see Fig. [TJ [N II] images). 

Although [N II] emission from the two jet-like 



features can be identified at slit positions 2 and 
5 (Fig. [6]), it can be better analyzed at slit posi- 
tions 7 and 8 which were specially selected to cover 
these structures. In the corresponding PV maps 
we can see the conspicuous emission arising from 
the two jet-like features. While the North jet-like 
feature emission is slightly redshifted, the South 
feature is slightly blueshifted. Both features show 
a small radial velocity, between and 5 km 
which suggest that they are moving almost in the 
plane of the sky. 

3. Discussion 

3.1. Morphokinematic structure and mod- 
eling 

As we have described, NGC 7354 shows four dif- 
ferent structures: two large scale structures, the 
outer and the inner shells, the last having two po- 
lar caps located on its symmetry axis; a number of 
low-ionization bright knots lying about the equa- 
torial plane, and two jet-like features located close 
to the North and South of the nebula, slightly in- 
clined respect to the inner shell symmetry axis. 
Our high dispersion spectra indicate that (a) the 
inner shell is expanding with a velocity which is 
smaller than that of the outer shell, (b) most of 
the bright equatorial knots are moving at veloci- 
ties closer to the outer shell, and (c) the jets' pro- 
jected velocity is very small. 

In order to test the overall morphology and 
kinematic structure of the nebula we have used the 
interactive three- dimensional (3D) mode ling tool 



SHAPE Ver. 2.0 (ISteffen fc Lopez I l2006l ) . SHAPE 



produces synthetic images and PV diagrams which 
can be compared directly with our observed CCD 
images and PV maps. We have considered two 
ellipsoidal structures (inner and outer shells), 
two spherical sections (polar caps), several small 
spheres with different diameter (bright knots) and 
two geometrically thin cylinders of 2" in diameter 
(jet-like features). All these geometrical struc- 
tures were slightly modified (smoothered, length- 
ened and/or rounded) in order to match the overall 
looking of our [N II] images. Besides, all the struc- 
tures were placed on the proper position to get 
the corresponding observed radial velocity from 
our PV maps. Fig. [7] shows all the geometrical 
components used to construct our shape model. 
Each one of the main structures was analized 
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separately in order to obtain its spatial velocity 
independently. Due to projection effects we are 
unable to distinguish neither if the low-ionization 
knots are located between the inner and outer shell 
nor if they are lying on the foreground or on the 
background side of the nebula. Thus, we have as- 
sumed that the bright knots are located on the 
outer shell wall, according to their individual ve- 
locity shown in our PV map. It is important to 
remark that the inclination angle of the different 
structures is unknown and SHAPE does not solve 
it. Final geometrical and kinematic parameters 
are shown in Table El Some of the inner shell ba- 
sic param^terscanbe compared to those found 



bv lHaiian et al. 



()2007() with an extended prolate 
ellipsoidal shell model. While our derived inclina- 
tion and position a ngles are cons i stent with the 
ones determined bv lHajian et al. ( 2007 ). our ex- 
pansion velocity value is larger than theirs. How- 
ever, as they have found in all the cases exam- 
ined in their study, the [N II] gas has a larger 
expansion velocity than the [O III] gas. Then, 
since we have determined the expansion velocity 
through the [N II] line emission, a difference in 
expansion velocities is expected. Distance esti- 
mations found in the literature lie around 1.2 kpc 



(ICahn et al.lll992l: IZhang Ill995l : IPhillips 12004 be- 
ing 1.5 kpc the largest estimation (jSabbadin 198^) 
and 0.88 kpc the smallest one (jPauh .19821 . We 
have adopted the largest value of 1.5 kpc in order 
to derive an upper limit for the kinematical ages. 
However, using a distance value of 1.2 kpc implies 
a decrease in the kinematical age estimation of 
only 20%. But most important, distance uncer- 
tainty does not modify the proposed chronological 
sequence of structure formation in our model. 

From this model, we have obtained synthetic 
images and PV diagrams which strongly resem- 
ble the observed ones, see Fig. [8] for an example 
of this good match. Thus, from the proposed 3D 
model we have not only reproduced qualitatively 
the overall morphology of NGC 7354 but we have 
also quantitatively reproduced the kinematic be- 
haviour of each structure. 

3.2. Formation of NGC 7354 

At present, binary star models have been suc- 
cesful in explaining the origin and process re- 
sponsible for the equatorial density enhacement 
required in the interacting stellar wind (ISW) 



model, or rather the generalized ISW (GISW) 
model, to explain early- and middle-type elliptical 
and even bip olar large scale structures observed 



PNe fe.g. iBalick et al.l Il987l: iMellema I Il995t 



Mellema fc Frank 1 119951 ). On one hand, common 
envelope (CE) evolution in close binaries has been 
proposed to be responsible for the expanding slow 



wind torus needed in GISW 



Rasio fc 



1995 iFrank et alJ Il99d : iTerman fc Taam I Il996t 
Sandauist et al. Ill998l) . On the other hand, ac 



ilVlO 



cretion disks in binary systems seem to account for 
the narrow waist bipolar morphologies and even 
the jets present in PNe. Details may vary on each 
binary model in order to explain individual mor- 
phologies and structures, but in all of them the 
main idea is that of a j et launched by the central 
star, or its companion ( Soker fc Rappaport 2000l: 



Garcia- Arredondo fc Frank 2004; Sahai et al. 



2007t iDennis et al.l l2008l lAkashi fc Sokeij l2008l 
among others). 

Since at first glance NGC 7354 looks like a 
rather elliptical nebula, one may think that the 
single star models (ISW or GISW) can explain the 
presence of the outer and inner shells, and even the 
formation of the two jet-like feature s observed in 



our [N II] image (jBalick et al.lll987^ . However, it 



is unable to explain why all the different structures 
show different position angles (PA's) on the sky. 
This characteristic may indicate that the direc- 
tion of ejection varies with time, and more likely 
that all the features were formed as independent 
events. Although nowadays there is no evidence 
of NGC 7354 possessing a binary nucleus, one may 
turn around and take a look at binary star mod- 
els to try to understand in a consistent way the 
observed morphology of NGC 7354. 

It has been proposed that the GISW model cou- 
pled with the predictions of the CE evolution the- 
ory can account for the formation of elliptical and 



even bipolar PNe morphologies (jSoker fc Livio 
19941 ). Then, we could succesfully explain the 
outer and inner shells in NGC 7354 based on 
this mixed model in the following way. After 
the CE phase has ended, i.e. the envelope has 
been ejected, the binary may become closer, mass 
transfer take place, an accretion disk is formed 



and a jet may be launched (jSoker fc Livio 1119941) 



If we assume that the subsequent evolution is sim- 
i lar to those described in t he models for pPNe 
(|Soker fc Rappaport I I2OOOI : iLee fc SahaTI l2003t 
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Dennis et alj booa lAkashi fc Sokeil 120081 ) it is 
possible that a jet arise from an accretion disk. 
Then, the accretion disk theory would be able 
to explain also the two jet-like features observed 
in NGC 7354. Moreover, some of these models 
even propose the formation o f an expanding dense 



ring in the equatorial plane (jSoker fc Rappaport 



2000) which in the case of NGC 7354 might be 
related to the equatorial bright knots. As the 
nebula evolve, a combination of the processes and 
mechanisms just mentioned might take place. 

With this in mind, we suggest a possible quali- 
tative scenario for the formation of the main struc- 
tures present in NGC 7354. A slow wind is lost by 
a central binary system ongoing a CE phase, pro- 
ducing the mild elliptical outer shell. When the 
AGB fast wind begins, it finds an equatorial den- 
sity enhaced enviroment forming the bright ellip- 
tical rim or inner shell. At this point the CE has 
been completely ejected and mass transfer from 
the secondary (main-sequence star) into the pri- 
mary (white dwarf) will form an accretion disk 
and eventually a couple of jets will be launched. 
Once the jets have appeared, they will "push" 
the polar ends of the inner shell, detaching them 
from it. Eventually, the two jets would break the 
caps and escape from the main body of the neb- 
ula. Along with the processes mentioned, the bi- 
nary system may be precessing and each of the 
structures would show a different PA on the sky. 
In the case of the two jets, precession may cause 
that they impinge on the polar caps with a cer- 
tain angle, i.e. not along the symmetry axis of the 
caps, breaking them apparently on the "base" of 
them. According to our derived kinematical ages, 
the chronological sequence of formation seem to be 
consistent with the above description. The outer 
shell is the oldest structure present in the nebula 
with an age of about 2500 yr. The next younger 
structure would be the inner elliptical shell with 
~1600yr. Kinematical ages for the two jet-like 
features indicate that they are coeval and both 
were formed almost at the same time as the inner 
elliptical shell. However, we should notice that 
the age was derived assuming that they have been 
moving with the same velocity (60 kms^^) since 
their launch. According with the theory of forma- 
tion of jets this low initial velocity would be very 
unlikely, since they arise from accelerated flowing 
material. Then, more likely the jets possessed a 



higher initial velocity but they have been restrain 
along its way. Maybe the main cause of this dec- 
celeration was the interaction of the jets with the 
polar caps. 

Regarding the equatorial low-excitation bright 
knots, as we have mentioned above they may be 
related to the dense expanding ring formed around 
the binary system. The bright knots show veloci- 
ties ranging between the inner and outer shell ex- 
pansion velocities. This may indicate that they 
are moving with a similar velocity as the par- 
ticular ambient gas in which they are embed- 
ded. We may compare these bright knots with 
the microstructures observed in NGC 2392 (es- 



kimo nebula, O'Deh et al. 1120021) and NGC 7662 



(jPerinotto et al. I l2004bl ) . Although each one of 



these nebulae is seen with a different view an- 
gle (NGC 2392 is seen pole-on and NGC 7662 is 
seen edge-on) both of them show various low- 
ionization knots located beyond the inner struc- 
ture and mainly on the outer envelope. If we com- 
pare the general morphology of NGC 7354 with 
that observed in NGC 2392 and NGC 7662 one 
can see that the three of them show an ellipsoidal 
inner structure surrounded by on outer envelope 
with low-ionization bright knots located on the 
outer structure plus jet-like features, interpreted 
as FLIERs in the case of NGC 7662. In the case 
of NGC 7354 since it is observed with an edge-on 
view angle and it is inclined with respect to the 
plane of the sky, we are not able to observe the 
real morphology of the equatorial brights knots. 
However, one can imagine that if one could get 
a pole-on view, i.e. with the inner shell major 
axis aligned with our line of sight, one would see 
bright knots located in a zone between the inner 
and outer shell very alike the arragement of knots 
with tails observed in the eskimo nebula. There- 
fore, having all these elements together, we suggest 
that, as in the cases of NGC 7662 and NGC 2392, 
the equatorial bright knots observed in NGC 7354 
are not spherical structures but knots with tails, 
or even filamentary structures. This last sugges- 
tion is strongly supported by the high resolution 
[N II] HST image. Fig. El 

A full model with MHD simulations is beyond 
the scope of the present study, but we are aware 
that it could clarify some aspects of our work. 



8 



4. Conclusions 

We have carried out a detailed morphological 
and kinematical analysis of the planetary nebula 
NGC 7354. In addition, we have derived the phys- 
ical conditions and elemental abundances in 21 re- 
gions of the nebula. The main conclusionas of this 
work can be sumarized as follows. 

Physical parameters of all the different struc- 
tures show that there are slight variations, both in 
electron density and electron temperature, within 
the nebula. Considering our error estimates, while 
density seems to slightly increase inward from the 
outer shell border to the interior of the inner el- 
liptical shell, temperature values seem to slightly 
increase in the same direction. Equatorial bright 
knots are clearly denser that both shells with tem- 
peratures in the typical range for ionized gas. 
Both jet-like features present a quite low density 
value and temperature was only determined for 
the South jet-like feature. All our derived physi- 
cal parameter values are consistent with the ones 
expected for radiatively excited gas. Local extinc- 
tion values within the nebula show a slight increas- 
ing gradient going from South to North of the neb- 
ula and from West to East along the location of 
the equatorial knots. 

Based on our kinematical data we have ob- 
tained a model that consistently reproduce the 
overall as well as the detailed morphology and 
PV maps of the nebula, using the 3D interactive 
tool SHAPE. The final geometrical components in- 
cluded in the model are: i) two shells with simi- 
lar inclination angle (respect to our line of sight) 
but different position angle of the projected semi- 
major axis; ii) two semispherical caps placed on 
the top of the inner elliptical shell; iii) thirteen 
spherical knots with different sizes, whose veloc- 
ities spanned between the corresponding velocity 
values for the outer and inner shells, placed at dif- 
ferent positions on the outer shell surface accord- 
ing to the optical image and observed spectra; and 
iv) two geometrically thin cylinders corresponding 
to North and South jet-like features. 

Finally, although there is no evidence of bina- 
rity in NGC 7354 at present, we suggest a quali- 
tative scenario for the formation of the different 
structures in the nebula based on the theory of 
CE evolution and formation of accretion disks in 
binary systems found in the literature. 
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Fig. 1. — Grey-scale representations of the images of NGC 7354 in the hght of Ha (left panels) and [N II]A6584 
(right panels). Two different contrasts have been used to show up the faint and bright nebular structures. 
Unsharped masking images are also present (lower row) in order to see details of the nebula. 
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Fig. 2.— Grcy-scalc representation of the F685N HST/WFPC2 image of NGC7354. The same features 
described from ground-based images can be identified here. The lower frame shows a lower contrast view of 
the equatorial region of NGC 7354. 
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Fig. 3. — Grey-scale and contour uniform-weighted map of NGC 7354 at 3.6 cm radio continuum. Contour 

levels are 8, 18, 35, 60, 100, 150, 200, 270, 350, and 400 times 144 ^Jy beam-\ the rms noise in the map. 
The half-power beam width (8'.'0 x 6'.'5, PA —86°) is shown in the upper right corner. 
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Relative RA (arcsec) 

Fig. 4. Slit positions (A,B,C,D) used for low resolution spectroscopy superimposed on the unsharped 
masking [N Il]A6584 image. Along each slit position several regions (1,2,3,.-.) were selected to analyze 
specific structures of the nebula. 
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Fig. 5. Slit positions used for high dispersion spectroscopy superimposed on the nnsharped masking [Nil] 
image. The slits are numbered from 1 to 8. Slit numbers 7 and 8 are oriented at PAs 25°and 13°, respectively. 
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Fig. 6. Grey-scale position- velocity maps of the Ha (upper panels), [N II]A6584 (middle panels) and 
He II A6560 (lower panels) emission lines derived from the long-slit spectra at slit positions 1 to 8 (indicated 
in the Ha panels for each column (see Fig. 5). The vertical and horizontal dashed lines in each panels indicate 
the systemic velocity of the nebula and the declination of the central star, respectively. North is up in each 
panel. 
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Fig. 7. — Upper panels show two sets of geometrical components used in the shape model. Lower panels 
show all components seen from two different angles, face-on (lower left panel) and pole-on (lower right panel). 
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Fig. 8. — Comparison of observed data and the 3D model obtained from SHAPE. Upper panels show the 
observed [N II] image (left panel) and the shape image (right panel). Lower panels show the observed PV 
map of the [N Il]A6584 line at slit position 7 (left panel) and the corresponding SHAPE synthetic PV map 
(right panel). 
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Table 1 

Dereddened emission-line fluxes (normalized to HB = 100). 
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Table 1 
— Continued 
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Ho I 


5876 


-0.215 


16.2 


16.7 


12.2 


10.5 


10.3 


8.8 


14.1 


16.2 


16.0 


15.8 


[O I] 


6300 


-0.282 


0.4 


1.6 


2.6 


1.8 


2.4 


0.8 


0.3 


0.4 




1.0 


[S HI] 


6312 


-0.283 


2.1 


1.8 


2.2 


2.3 


2.3 


2.0 


1.9 


1.5 


2.2 


1.5 


[OI] 


6364 


-0.291 






0.9 


1.0 


0.7 


0.2 










[N II] 


6548 


-0.318 


8.0 


20.8 


28.2 


22.1 


26.8 


10.5 


7.0 


18.4 


5.8 


8.8 


Ha 


6563 


-0.320 


280.2 


284.2 


283.5 


283.5 


283.9 


282.0 


283.5 


282.7 


285.0 


283.0 


[N II] 


6584 


-0.323 


24.5 


65.2 


82.8 


66.6 


82.8 


31.8 


21.0 


58.4 


19.5 


28.2 


He I 


6678 


-0.336 


4.1 


4.4 


3.5 


3.5 


3.2 


2.7 


4.2 


4.6 


4.5 


5.1 


[S II] 


6717 


-0.342 


1.9 


3.5 


3.6 


2.8 


3.4 


1.8 


2.2 


4.2 


1.9 


2.1 


[S II] 


6731 


-0.344 


2.5 


5.6 


5.0 


4.1 


5.1 


2.7 


2.2 


5.3 


2.5 


3.4 


Ho I 


7065 


-0.387 


3.6 


4.2 


3.5 


3.3 


3.3 


2.6 


3.3 


4.2 


4.5 


3.7 


[Ar III] 


7135 


-0.396 


16.0 


19.6 


18.0 


17.3 


17.2 


15.2 


14.5 


18.7 


15.9 


15.8 


[Ar IV] 


7236 


-0.409 


0.5 




0.9 


1.1 


1.0 


1.3 










Ho I 


7281 


-0.414 


0.7 




0.6 


0.5 


0.4 


0.6 






1.4 


0.4 


[O II] 


7320 


-0.419 


1.0 


1.6 


2.3 


2.0 


1.9 


1.2 


0.5 


2.9 


1.1 


0.6 


[O II] 


7330 


-0.420 




1.4 


2.1 


2.1 


2.3 


1.3 


1.4 








c(H/3) 






2.01 


2.20 


2.41 


2.24 


1.95 


2.02 


1.71 


1.78 


1.83 


1.83 


log /(H^) 






-15.801 


-16.357 


-16.056 


-16.224 


-15.203 


-15.313 


-15.399 


-15.967 


-15.874 


-15.913 



20 



Table 2 
Derived physical parameters. 



Region 




[K] 




[cm-3] 






[Nil] 


[O III] 


[S II] 


[CI III] 


[Ar IV] 


A 1 






1 1 zl 1 4-1 flSfl 






A O 


ioyuy±:il4y 




ZU4oitoOU 


0D40dl ZDDU 


i^i=;f^Q4- "1 Q T o 

ooooit loiy 


A3 




1 2281 +2(S2 


262*^+281 


402'H4- 1888 




A4 






1982±5o7 


1033^952 


8812-'-1201 


A ''i 


1 - 1 o",-|- 1 (ion 


12135i i 17 


1 iO ii ioS 




2;-; 11 — 15(S7 


Bl 




13130±700 


751±267 






B2 




14346±519 


1799±610 






B3 


17690±1468 


15523±356 


2402±515 






B4 




13885±263 


3106±555 


7708±4389 


1062ib767 


B5 




13646±286 


2216±415 


15086±9695 




B6 


12833±567 


13005±298 


2700±356 




2986ibl423 


B7 




14334±956 


1910±636 






CI 


10724±2902 




3451±4044 






C2 


11312±785 


10363d=1162 


2130±629 






C3 


11329±1438 




2354±1194 






C4 


11006±481 


12478±755 


2715±551 






C5 


12662itl271 


15454±1039 


3311±1290 






C6 




11298±1564 


664±910 






Dl 




12066±984 


1603±575 






D2 






1574±725 






D3 




11825it835 


4102±2481 







Table 3 
Derived ionic abundances. 



Region 












Ion 














He+ 


He++ 




o" 




Q + 2 


S + 




Ar+2 












[IxlO^'^j 


[lxlO-'5] 


[IxlO-'^j 


[IxlO-"] 


1x10"^] 


[IxlO-'^j 


[IXlO-'^j 


[IxlO-**] 


[1x10"''] 


Al 


0.117 


0.003 


1.38 


10.12 


2.98 


4.89 


2.83 




1.30 






A2 


0.073 


0.046 


0.10 




0.18 


1.10 


0.37 


0.80 


0.55 


0.49 


0.69 


A3 


0.074 


0.042 


0.60 


0.74 


1.11 


2.69 


1.63 


3.17 


1.15 


0.91 


1.00 


A4 


0.058 


0.062 


0.13 


0.15 


0.49 


1.90 


0.42 


1.47 


0.82 


1.04 


1.06 


A5 


0.099 


0.020 


0.19 


0.50 


0.28 


1.40 


0.47 


0.81 


0.65 


0.46 




Bl 


0.112 


0.007 


0.42 




3.98 


4.76 


1.45 


6.68 


1.59 






B2 


0.082 


0.038 


0.26 


0.72 


1.60 


5.15 


0.87 


5.54 


1.54 


1.32 




B3 


0.062 


0.048 


0.08 


0.11 


0.13 


0.94 


0.30 


0.74 


0.48 


0.37 


0.39 


B4 


0.069 


0.055 


0.25 


0.35 


1.86 


4.14 


0.99 


3.99 


1.39 


1.54 


2.50 


B5 


0.061 


0.061 


0.23 




1.71 


3.98 


0.82 


4.88 


1.29 


1.60 


2.90 


B6 


0.074 


0.045 


0.46 


0.86 


0.57 


2.13 


1.12 


1.89 


0.84 


0.67 


0.85 


B7 


0.110 


0.007 


0.48 


0.74 


1.88 


4.84 


1.39 


4.89 


1.48 


1.03 




CI 


0.110 


0.001 


1.11 


2.46 


1.56 


3.29 


2.87 


3.08 


1.53 






C2 


0.083 


0.028 


1.21 


3.22 


2.13 


3.08 


2.06 


3.22 


1.26 


1.32 


0.88 


C3 


0.088 


0.032 


0.98 


3.41 


1.96 


3.14 


1.69 


3.27 


1.20 


1.64 


1.11 


C4 


0.080 


0.047 


1.31 


3.18 


2.46 


3.01 


2.35 


3.69 


1.27 


1.02 


1.07 


C5 


0.068 


0.051 


0.36 


0.63 


0.60 


1.89 


0.98 


1.92 


0.84 


0.81 


1.64 


C6 


0.101 


0.018 


0.41 


0.67 


1.34 


4.98 


1.15 


4.69 


1.33 


1.56 




Dl 


0.116 


0.005 


1.13 


0.69 


5.51 


4.88 


2.85 


3.51 


1.73 






D2 


n.11.3 


0,006 


0,38 




2,18 


4.65 


1.31 


5.21 


1.46 






L):! 


n.ll:i 


0.007 


0,,")7 




0,hT 


1.92 


2.10 


3. .11 


1. 1,-) 






Note.— 


-In regions A5, Bl and D3, He++ 


abundances shown are an upper limit. 
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Table 4 
Derived total abundances. 



Region 


TT 

Hc/H 


O/H 






N/H 




S/H 






Al / rl 


Al 


0.120±0.010 


8.72±0. 


,05 


8 


.39±0.35 


6, 


,67±1, 


,54 


6 


.38±0.10 


A2 


O.llSitO.OOl 


8.19±0. 


14 


7 


.94±0.28 


6, 


,41±0, 


14 


6 


.05±0.01 


A3 


0.116±0.001 


8.58±0. 


15 


8 


.31±0.21 


6, 


,88±0, 


,09 


6 


.35±0.01 


A4 


0.121±0.001 


8.50±0. 


,40 


7 


.94±0.60 


6, 


,63±0, 


,27 


6 


.30±0.01 


A5 


0.119±0.001 


8.21±0. 


,28 


8 


.02±0.59 


6, 


,36±0, 


28 


6 


.09±0.14 


Bl 


0.118±0.002 


8.73±0. 


11 


7 


.76±0.13 


7, 


,06±0, 


,05 


6 


.47±0.10 


B2 


0.120±0.002 


8.84±0. 


,08 


8 


.05±0.12 


7, 


,14±0, 


,04 


6 


.46±0.10 


B3 


O.llOitO.OOl 


8.15±0. 


18 


7 


.90±0.26 


6, 


,40±0, 


14 


5 


.95±0.01 


B4 


0.124±0.001 


8.81±0. 


,04 


7 


.94±0.07 


6, 


,97±0, 


,03 


6 


.51±0.01 


B5 


0.122±0.001 


8.82±0. 


,05 


7 


.94±0.07 


7, 


,07±0, 


,03 


6 


.51±0.01 


B6 


0.119±0.001 


8.48±0. 


15 


8 


.38±0.21 


6, 


,72±0, 


10 


6 


.21±0.01 


B7 


0.117±0.002 


8.72±0. 


16 


8 


.13±0.18 


7, 


,03±0, 


,07 


6 


.44±0.10 


CI 


0.111±0.006 


8.54±0. 


,89 


8 


.39±4.87 


6, 


,82±1, 


,90 


6 


.46±0.40 


C2 


0.111±0.003 


8.60±0. 


,30 


8 


.36±0.45 


6, 


,81±0, 


20 


6 


.45±0.03 


C3 


0.120±0.006 


8.61±0. 


,27 


8 


.31±0.74 


6, 


,83±0, 


,36 


6 


.49±0.04 


C4 


0.127±0.002 


8.65±0. 


,26 


8 


.37±0.33 


6, 


,86±0, 


14 


6 


.41±0.02 


C5 


0.120±0.003 


8.45±0. 


,55 


8 


.23±0.72 


6, 


,71±0, 


29 


6 


.27±0.02 


C6 


0.120±0.005 


8.76±0. 


,23 


8 


.25±0.97 


7, 


,07±0, 


,33 


6 


.40±0.10 


Dl 


0.122±0.007 


8.75±0. 


14 


8 


.06±0.16 


6, 


,77±0, 


,08 


6 


.51±0.10 


D2 


0.118±0.002 


8.70±0. 


,01 


7 


.94±0.13 


7, 


,03±0, 


,06 


6 


.44±0.10 


D3 


0.120±0.002 


8.72±0. 


13 


8 


.53±0.17 


7, 


,01±0, 


,09 


6 


.43±0.10 


PNo" 
TI PNo*" 
H 11"= 
Sun"* 


0.12±0.02 
0.13±0.04 
0.10±0.01 
0.09±0.01 


8.68±0.15 
8.65±0.15 
8.70±0.04 
8.66±0.05 


8.35±0.25 
8.72±0.15 
7.57±0.04 
7.78±0.06 


6.92±0.30 
6.91±0.30 
7.06±0.06 
7.14±0.05 


6.39±0.30 
6.42±0.30 
6.42±0.04 
6.18±0.08 



Note. — Except for He, all abundances relative to H are logarithmic 
values with l ogH — +12. Final rows are for comparison: (a) Aver- 
age for PNe dKingsb urghfc Barlow I |199^; (b) Avera ge for Type I PN e 
llKinssbureh & Barlow 199 1); (c) Aver age for H II regions dShaver et al. 1119831^ : 
the Sun (Grevesse, Asplun d Sauval I i2007i) . In regions A5, Bl, and D3, He 
abundances shown arc an upper limit. 



Table 5 

Geometrical and Kinematic Parameters from our SHAPE model. 



Geometrical Structure 


Size 

["] 


Expansion Velocity 
[km s"^] 


Inclination Angle i 
[dog] 


Position Angle tp 
[dcg] 


Kinematical Age 

[y] 


Outer Elliptical Shell 


16 X 15 


45 


94 


18 


2 500 


Inner Elliptical Shell 


10 X 8 


35 


94 


28 


1 600 


Scmispherical cap N 


14 


47 


94 


28 


1 000 


Somispherical cap S 


14 


47 


85 


203 


1 000 


Spherical Knots 


1 to 2 


25 to 44 








Thin Cylinder N 


7 


60 


95 


3 


1 900 


Thin Cylinder S 


13 


60 


85 


205 


1 900 



Note. — Size refers to major and minor axis in the case of elliptical structures, length for thin cylinders, and diameter in the 
case of the spherical knots and polar caps. All kinematical ages were determined assuming a distance to the nebula of 1.5 kpc 
jSabbadin I IT986 V 
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